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PRSGEDINO  PAOBi BLANK-NOT  FILMED 


I 


INTRODUCTION 


The  reader  is  assumed  to  be  familiar  with  the  axiByrnmetric  version  of 
the  SMITE  code  which  is  described  fully  in  BRL  CR-255.  It  is  also  assumed 
that  the  reader  is  familiar  with  the  Stanford  Research  Institute  (SRI) 
ductile  fracture  subroutine,  DFRACT,  and  bi'ittle  fracture  subroutine,  BFRACT. 
The  description  of  the  theory  of  fracture  used  in  these  routines  are  con¬ 
tained  in  AMMRC  CTR  75-2  and  BRL  CR-222.  . 

The  purpose  of  this  research  effort  was  twofold; 

(i)  First  to  modify  the  SMITE  code  so  that  the  ductile  and  brittle 
fracture  model  (that  of  SRI)  could  be  incorporated  into  SMITE. 

(ii)  Compare  results  obtained  with  the  modified  code  with  the  La- 
grangian  formulation  of  SRI. 

The  modifications  introduced  into  the  SMITE  code  to  incorporate  the 
fracture  subroutines  are  described  in  the  following  section.  Each  subroutine 
that  required  changes  is  described. 
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II.  FRACTURE  MODIFICATIONS  TO  SMITE 


DENSE  -  Since  pressure  is  now  carried  as  an  explicit  variable,  the 
pressure  on  the  free  surface  Is  set  to  that  necessary  to  produce  zero  normal 
stress.  The  Iteration  of  density  to  obtain  the  required  pressure  through 
the  solid  equation  of  state  Is  removed. 

FDIFF  -  Logic  is  added  to  call  IFRACT  and  hence  obtain  the  relative 
properties  from  the  fracture  subroutines  once  the  critical  pressure  has  been 
'Exceeded  at  a  point. 

IFRACT  -  This  is  the  routine  that  actually  incorporates  the  fracture 
subroutine  DFRACT  or  BFRACT.  This  is  accomplished  by  introducing  a  quasi- 
Lagranglan  calculation  Into  the  Eulerlan  code.  Once  the  velocities  at  the 
advanced  time  are  obtained  at  a  mesh  point  by  the  two  step  method,  they  are 
integrated  backwards  in  time  to  determine  where  that  mesh  point  came  from  in 
a  Lagranglan  sense.  Values  of  all  variables  at  the  initial  time  are  then 
assigned  to  this  point  by  linear  interpolation  from  surrounding  points.  The 
ductile  fracture  subroutine  DFRACT  or  the  brittle  fracture  subroutine  BFRACT 
may  then  be  called  to  update  these  values  to  the  original  mesh  point  at  the 
advanced  time. 

Values  for  strain  are  obtained  by  numerically  differentiating  the  veloci¬ 
ties  using  central  differences.  An  a  acceleration  at  a  point  with  index  i, 
for  instance,  may  be  obtained  by  differencing  velocities  at  points  i-1  and 
i+1.  If  the  Lagranglan  point  lies  between  mesh  lines  i  and  i+1,  acceleration 
at  i+1/2  may  also  be  obtained  by  differencing  velocities  at  points  i  and  i+1. 
Similar  differencing  is  used  to  obtain  B  accelerations.  Accelerations  at  the 
Lagranglan  point  are  then  evaluated  by  interpolation.  The  necessary  strains  are 
obtained  by  transforming  from  a,  6  to  z,  r  space. 


6 


INPACE  -  The  pressure  on  both  siSes  of  the  Interface  is  set  to  that 
necessary  to  produce  a  continuous  norm&l  stress.  The  iteration  of  density  to 
obtain  the  required  pressure  through  the  solid  equation  of  state  is  removed. 

INITAL  -  The  number  of  dependent  variables  is  changed  to  10  for  ductile 
fracture  and  22  for  brittle  fracture.  Seven  pareuneters  are  read  in  for  duc¬ 
tile  fracture  and  fourteen  for  brittle  fracture. 

INTRPL  -  Coding  is  added  to  reflect  the  additional  variables  about  the 
axis  of  symmetry. 

OUTPUT  -  Coding  is  added  to  print  the  additional  variables  at  points 
where  fracture  has  occurred. 

PLTOUT  -  The  output  of  points  that  are  in  the  plastic  region  is  re¬ 
moved  and  the  eimount  of  damage  to  fractured  points  is  substituted. 

PRS  -  The  solid  equation  of  state  is  changed  to  the  Mie-Gruneisen 
equation. 

P  =  (Cu  +  Du^  +  Su^)  (1  -  p-)  +  rpE 

where  C,  D,  S  are  constants,  T  is  the  Gruneisen  ratio,  E  is  the  internal 
energy,  P  is  the  density  and  u  =  p/p^-l.  In  the  next  section,  we  describe  the 
modifications  required  for  input  data  preparation. 
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III.  INPUT 


There  are  three  types  of  Input  data,  Integer  (I),  real  (R)  and  alpha¬ 
numeric  (A) ,  An  Integer  la  a  number  without  a  decimal  point  which  must  be 
right  justified  in  its  field.  A  real  number  has  a  decimal  point  which  may 
be  followed  by  an  exponent  of  the  form  Ej;N.  The  +N  represents  the  power  of 
10  by  which  the  number  is  to  be  multiplied.  The  +  may  be  omitted  if  N  is 
positive  and  the  E  need  not  appear  if  either  the  -f  or  -  is  present.  If  an 
exponent  is  present  it  must  be  right  justified  in  the  field.  Alphanumeric 
input  consists  of  exactly  the  punched  characters  and  blanks  appearing  in  the 
field.  A  description  of  the  necessary  input  data  for  SMITE  follows.  The 
card  number  refers  to  the  type  of  input.  If  more  than  one  card  is  necessary 
for  this  input,  there  may  be  several  cards  of  the  same  type  number. 
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The  set  of  cards  containing  the  input  parameters  for  the  individual  materials  follows 
This  entire  set  is  repeated  for  each  material  until  all  materials  have  been  defined. 


I  Ni  MaxUBia  n’xitMr  of  mesh  llnaa  allowed  in  axial  direction. 

6-10  I  NJ  HaxiBvxB  nxBber  of  nesh  lines  allowed  in  radial  direction 
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IV.  DUCTILE  FRACTURE  RESULTS 


A  tapered  flyer  plate,  Alumlnuin-1145,  impacting  on  a  flat  plate, 
Aluminuin-1145,  at  a  striking  velocity  V-.0163  cm/U-sec.  was  chosen  as  a  test 
c?‘se.  It  is,  except  for  the  cyla  drical  geometry,  similar  to  the  problem  run 
by  SRI  which,  in  that  case,  was  i  >ecified  in  plane  coordinates.  The  coordin¬ 
ate  mesh  chosen  for  the  computation  presented  here  is  not  the  finest  possible 
but  is  similar  to  the  SRI  mesh. 

The  target  had  a  mesh  spacing  Az«L/20  where  the  thickness  L=0.316  cm. 

The  spacing  of  the  mesh  in  the  radial  direction  Ar=h/40  where  the  height  of 
the  target  is  h=1.184  cm.  In  the  tapered  fly  r  plate  the  mesh  spacing  in 
z,  Az  =  0.158  cm./12)  the  radial  mesh  increment  is  Ar=0.592  cm./24. 

Tne  following  table  gives  the  parameters  (data  card  10b)  used  in  the 
nucleation  and  fracture  model: 

FRACTURE  PARAMETERS:  Aluminum-1145 


Growth  Constant 

-0.01 

Growth  Threshold 

-0.4  •  10^° 

Nucleation  Radius 

0.1  •  10^ 

• 

10 

% 

0.3  •  10 

Nucleation  Threshold 

-0.3  •  10^'^ 

"l 

-0.4  •  10^ 

The  results  of  the  computatio:.  for  the  value  of  the  parameters  given  in 
the  above  table  is  presented  in-  the  three  accompenying  figures.  Figure  1 
shows  the  geometrical  configuration  at  the  moment  of  impact,  t=0.  At  t=1.464 
(isec.,  shown  in  F.-'gure  2,  f'.d  >  j’-ows  i'.'s  presence  or  an  internal  frac¬ 

ture  centered  about  the  ce'  ter  line  of  the  plate.  This  fractured  region  re¬ 
presents  a  percent  fracture  volume  between  5  to  12  percent. 
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The  presenca  of  thla  fi'acture  island  Is  similar  to  that  obtained  by 
SRI,  however  their  calculation  was  carried  out  with  the  brittle  model  for 
a  flyer  plate/target  configuration  with  the  material  properties  of  Armco 
iron.  As  a  result,  due  to  the  higher  growth  rates  of  this  brittle  material, 
the  fracture  Island  was  of  much  greater  extent,  a  reasonable  result  (see  the 
Table  below).  The  last  figure  shows,  at  ta2. 198V -sec. ,  that  there  is  vir¬ 
tually  no  continued  growth  of  the  fracture  Island. 

FRACTURE  PARAMETERS!  Armco  Iron 


Growth  Constant 

-0.12  .  10 

Growth  Threshold 

0. 

-2 

Nucleation  Radius 

0.4  .  10 

NO 

0.5  .  10® 

Nucleation  Threshold 

1.12  .  10^^ 

-7.4  .  10® 

The  brittle  fracture  routine  was  substituted  for  the  ductile  fracture 
routine,  and  the  axisymmetrlc  flyer  plate  problem  was  run  with  Armco  iron 
properties  rather  than  Aluminum-1145  properties.  At  this  point,  we  have  very 
limited  results.  The  main  difficulty  is  a  convergence  problem  in  the  brittle 
fracture  subroutine.  We  have  been  conversing  with  the  scientific  staff  at 
SRI  who  are  responsible  for  the  brittle  fracture  model.  As  a  result  of  their 
suggestions,  there  has  been  improvement  in  the  convergence  of  the  iterative 
procedure  in  BFRACT  but  not  a  coi!?>lete  solution  to  the  problem. 
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SUMMARY  AND  RECOMMENDATIONS 


The  feasibility  of  adapting  the  SRI  Lagrangian  fracture  subroutines  for 
inclusion  in  the  SMITE  Eulerian  code  has  been  demonstrated.  The  ductile 
fracture  subvoutine  was  successfully  run  with  reasonable  results.  The 
brittle  fracture  subroutine  was  run,  but  aborted  due  to  convergence  problems., 
SRI  was  consulted  about  the  problem  and  recommended  corrections  were  received. 
However,  the  corrections  were  not  received  until  after  the  funds  under  the 
present  contract  were  expended.  Those  corrections,  therefore,  have  not  been 
included  in  the  code. 

It  is  recommended  that  the  following  additional  items  be  considered  by 

BRL: 

1)  Inclusion  and  testing  of  the  corrections  recommended  by  SRI  with  addi¬ 
tional  consultation,  if  necessary. 

2)  More  extensive  testing  of  both  the  ductile  and  brittle  fracture  sub¬ 
routines  to  determine  if  any  further  corrections  are  necessary. 

3)  Further  modificati<^n  of  printed  and  graphic  output  to  present  re¬ 
sults  in  a  form  suitable  to  the  requirements  of  BRL. 
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